
SYSTEMS AND METHODS OF DETERMBSflNG MOTION TOOLp ^ 
PARAMETERS IN BOREHOLE LOGGING ^ ^ 



LP OF THE INVENTION ^. S 

The present invention generally relates to borehole logging ai^moreg 
particularly, to systems and methods for determining parameters of the motion of a logging 
tool in a borehole and compensating for this motion in the logging data. 

10 BACKGROUND OF THE INVENTION 

Knowledge of the motion parameters of a logging tool relative to a borehole 
is important for obtaining accurate downhole measurement signals. Thus, for example, a 
shift in the sensitive volume of the logging tool caused by lateral motion of the tool relative 
to the borehole can distort the received signal and produce inaccurate measurements. Such 
15 distortions can be especially problematic in logging-while-drilling (LWD) and 

measurement-while-drilling (MWD) environments, where the tool itself is subjected to 
severe vibration. In some cases, the measurements may have to be completely discarded, as 
for example when a stick-slip condition occurs (where the drill bit stops rotating 
momentarily because of high friction and then rapidly accelerates when the moment exerted 
20 by the drill pipe exceeds the static friction). Clearly, it would be beneficial if the tool 
operator had access to information about the motions of the tool, so that measurements 
made during strong lateral and axial motions are discarded, or not even initiated. 

Additionally, in many cases, it is important to select optimal activation times 
for the logging measurement and, if possible, to enable correction of the received signal 
25 based on the motion parameters of the logging tool. In such applications, it is necessary to 
accurately determine the lateral tool velocity of the tool in real time. 

In the simplest system for measuring the lateral tool velocity relative to a ^ 
borehole, two mutually orthogonal accelerometers can be placed on the tool axis, such that 
they are sensitive to the lateral acceleration. However, such placement is generally not 
30 possible in downhole tools because of design constraints, in particular owing to the need to 
have an open space within the center of the tool for a mud flow channel. 

As such, in prior art systems for determining the lateral velocity of a drilling - 
tool, two pairs of accelerometers are attached to diametrically opposite sides of the tool, 
usually near the surface of the tool. See, for example, co-pending appUcation Ser. No. 
35 i ^/di^ f2, Jiiud J miUdi l y ' 1 J, IW b K assigned to the assignee of the present appUcation. This 
t^ .S. P. application is incorporated herein by reference for all purposes. The accelerometers 
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together provide radial acceleration components, a^j and a,2, and tangential acceleration 
components, a^j and a^^ of the tool. Since the accelerometers rotate with the tool, their 
measurements are in the reference frame of the rotating tool, i.e., the rotating frame. Given 
their opposite placement, the accelerometer pairs register eqiial but opposite accelerations 

5 due to lateral tool motion and equal radial (centrifugal) as well as angular accelerations due 
to tool rotation. The radial arid tangential forces due to tool rotation are compensated for 
the opposite accelerations by subtracting the reading of one accelerometer from the reading 
of the diametrically opposite one (%2 is subtracted from a^j and is subtracted from 
The remaining signal is twice the actual lateral tool acceleration in the directions of a^i and 

10 a^i, respectively, as seen in the rotating frame. The acceleration components compensated 
for the centrifiigal and angular accelerations are therefore given by the expressions: 

a, = (a^i - ay2y2, for the radial tool acceleration; and 
a^ = (a^i - at2)/2, for the tangential tool acceleration. 

15 

The lateral velocity is calculated by integrating the above acceleration 
components. There are two main problems associated with this prior art approach. First, 
the signal measured by the accelerometers will also contain a gravitational component if the 
tool orientation is not vertical. The magnitude of the gravitational component is Gsina, 

20 where a is the angle of tool inclination relative to vertical and G is the gravitational 

acceleration constant. The frequency of the gravitational component is related to the tool 
rotation. Tool tilt away from vertical is sensed by the accelerometers and, thus, introduces 
an inaccuracy in the lateral tool acceleration readings. 

Commonly, the gravitational acceleration component is removed from the 

25 signal by employing a high pass filter. The filter cut-off frequency is set to separate 

frequencies of the gravitational component modulated by the tool's rotation from the higher 
frequencies assimied to be caused by the tool's lateral motion. This technique, however, is 
not effective if the tool's rotational rate is high or not constant, for example, in a stick-slip 
situation, gravitational acceleration components are generated within the band of those 

30 related to the tool's lateral motion. 

The second problem occurs because the accelerometers, which are placed on 
the tool, measure the tool's lateral velocity in the tool frame of reference rather than the 
desired borehole frame of reference. 
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With reference to Fig. IB, while the motion parameters a^ and are provided 
in the rotating reference frame, it is desirable to determine the corresponding motion 
parameters of the tool in the borehole reference frame XYZ. It will be appreciated that if 
the tool does not rotate, then the tool-reference parameters a^ and a^ are equivalent to the 

5 borehole reference parameters a^ and ay, and no conversion is necessar3(. If, however, the 
tool rotates, then and a, are different than a^ and ay, and conversion to the borehole frame 
of reference will be required. Similarly, it will be necessary to convert velocity components 
Vj and v^ to v^ and Vy , corresponding to the borehole reference frame. 

Obtaining accurate lateral tool velocity is important to ensure that the 

1 0 accuracy of NMR porosity measurements does not degrade by more than about 5%. For 
example, the lateral displacement of the tool's center of gravity should be limited to about 
0. 1 mm relative to the borehole within a measuring time frame of 500 (^seconds. In 
practice, it is desirable that the lateral tool velocity should not exceed 0.2 m/sec during a 
typical NMR reading. Tool displacements greater than about 0.25 mm may introduce a 

1 5 system error associated with phase shift of the NMR echo. In addition to systemic error, the 
signal-to-noise may also degrade. 

Therefore, there is a need to provide a system and method for accurately 
determining the lateral tool velocity and overconaiing the deficiencies associated with the 
prior art. By knowing the tool's velocity, the NMR signal may be corrected. Additionally, 

20 along with velocity information, an xmcertainty estimator can be calculated to provide 
confidence levels of the measurements. 

SUMMARY OF THE INVENTION 

The present invention is directed to a system and method for providing 
25 lateral tool velocity measurements corrected for a gravitational acceleration component, and 

for lateral velocity measurements in the borehole frame of reference. The invention allows 

correction of either one independently. 

The system of the present invention comprises two pairs of accelerometers, a 

pair placed on opposite radial sides of the tool's rotational axis. Each of the two pairs of 
30 accelerometers detects both radial and tangential tool acceleration components. The system 

fiirther comprises: two magnetometers placed orthogonally with respect to each other to 

detect the tool's magnetic phase (also called magnetic tool face); means for sampling the 

accelerometer and magnetometer signals; and means for real time processing and 

calculation of the corrected lateral tool velocity. 

35 
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The method of the present invention includes: concurrently obtaining the 
tool's instantaneous radial and tangential accelerometer data as well as the tool's 
instantaneous magnetic phase data; determining the tool's gravitational phase data, and 
using this data to provide a corrected lateral tool velocity. 

5 Accordingly, in one aspect, the invention is an apparatus for making 

borehole measurements using a logging tool, comprising: at least two accelerometers 
measuring motion parameters of the tool in a borehole, the motion parameters being 
measured in a tool reference frame; two or more magnetic induction sensors meastuing 
magnetic induction dependent on the orientation of said accelerometers in relation to the 

10 gravitational field; and a processor computing the motion parameters of the tool in a 
borehole reference frame by taking into accoxmt the orientation of said accelerometers in 
relation to the gravitational field. 

In another aspect, the invention is a method for making borehole 
measiirements with a logging tool moving along a borehole subjected to a gravitational 

15 field, comprising: (a) measiuing motion parameters of the tool in a tool reference frame; 
(b) estimating components of the measured motion parameters due to the gravitational field 
in the borehole; (c) computing a set of corrected motion parameters, the step of computing 
comprising removing the estimated components due to the gravitational field; and (d) 
providing at least one measurement along the borehole with the logging tool, said at least 

20 one measurement taking into accoimt the computed set of corrected motion parameters. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention is best understood with reference to the attached 
drawings, in which: 

25 Figure lA shows a schematic diagram of the processing system; 

Figure IB depicts a generalized view of the rotating drill placed within a 

borehole; 

Figure 2 A and 2B show cross-sectional views of a rotating dbrill tool with 
accelerometers in alternative placements; 
30 Figure 3 shows the variation of the magnetic tool phase obtained in 

accordance with the invention with two magnetometers placed on the tool; 

Figure 4 is a block diagram showing the processing required to detmnine the 
gravitational component of the acceleration signals in a preferred embodiment; 

Figure 5 is a schematic diagram illustrating the processing steps for 
35 calculating velocity using a multiple window approach; 
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Figure 6 is a graph of instantaneous angular velocity with low vibration; 
Figure 7 is a graph of instantaneous angular velocity in stick-slip conditions; 
Figure 8 is a graph of the spectral power density of lateral acceleration under 
stick-slip conditions; 

5 Figure 9 is a graph of a computer synthesized frequency modulated signal in 

a simulated stick-slip condition; ' 

Figure 10 is a graph of the spectrum of a computer synthesized frequency 
modulated signal; 

Figure 11 is a graph of a Fourier transform of the accelerometer signal during 

10 bit whirl; 

Figure 12 is a graph of simulated accelerations diuing collision of the drill 
collar with a rigid barrier; 

Figure 13 is a graph of peak acceleration versus initial or impact velocity; 
Figure 14 is a graph of the effect of limited bandwidth on the acceleration 
15 signal after it passes through a 500 Hz low pass filter; 

Figure 15 is a graph of velocity error due to group delay of a low-pass filter; 
Figure 16 is a graph depicting velocity error due to group delay vs. 

bandwidth; 

Figure 17 is a diagrammatic view of the tool axis in the borehole, the tool 
20 having a drill collar and stabilizer; 

Figure 18 illustrates NMR antenna magnet influence on magnetic sensor 

measurements; 

Figure 19 is a graph of the relationship between the magnet-magnetometer 
distance and the By component; 
25 Figure 20 is a graph showing the sensitivity of the By magnetometer 

component versus deformation of the magnet-magnetometer assembly expressed in magnet- 
magnetometer distance; and 

Figure 21 is a graph of velocity uncertainty combining all three 
magnetometer-related effects. 



35 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The structure and function of the preferred embodiment can best be 
understood by reference to the attached drawings. Where the same reference numerals 
appear in muhiple figures, the nimierals refer to the same or corresponding structure in 
S those figures. 

♦ . 

A, The System of the Invention 

In a preferred embodiment, the system of the present invention for 
determining the lateral velocity of a drilling tool comprises at least two pairs of 

10 accelerometers placed opposite each other across the tool axis of rotation, a first and second 
magnetometer placed on the tool to provide detection of the tool's rotational magnetic phase 
relative to the earth's gravitational field, an interface for reading signals firom the 
accelerometers and magnetometers, and a data processor for providing a corrected value of 
lateral tool velocity in the borehole firame of reference. 

15 As shown in Figure 1 A, in a preferred embodiment of the system, the 

interface for reading data signals comprises electrical or electronic circuitry (not shown) that 
is associated with the use of accelerometers and magnetometers, as known by those skilled 
in the art. The data processor generally is implemented using a microprocessor 1 1 which 
runs a software program 10 for accepting digital data signals firom the accelerometers and 

20 magnetometer 13 and processing the raw data signals for instantaneous conversion. The 
interface includes an analog to digital (A/D) signal converter 9 used to convert raw data 
signals into digital signals. Appropriate software programming 10 is used to process and 
calculate corrected lateral tool velocities using the data. The system fiirther provides an 
NMR output trigger 14 and a terminal output IS to a monitor, printer, or similar output 

25 device. An optional data storage 12 may also be used for storing raw or processed data for 
later retrieval and processing. 

Examples of equipment that can be used in accordance with the present 
invention include the digital signal processor model ADSP21060 SHARC chip firom Analog 
Devices, which possesses a six chaimel A/D converter to accommodate acceleration and 

30 magnetic signals, micro-machined siUcone accelerometers firom Silicon Designs with a 
bandwidth of SkHz and niagnetoresistive sensors (magnetometers) from Honeywell (Model 
HMC1002). Also included in a preferred embodiment of the system are conmiunication 
channels (not shown) to control the digital processor and to retrieve real time velocity 
information, and a digital synchronization line to trigger the NMR measurements when the 

35 instantaneous lateral tool velocity is low. 
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Referring to Fig, IB and Figures 2 A and 2B, in accordance with tiie present 
invention, the XYZ reference frame is related to the borehole 1, where the tool is placed. 
The (r-t-a) reference frame, on the other hand, is related to tool 2, which in general may be 
rotating. Accordingly, logging-while-drilling (LWD) and measurement-while-drilling 

5 (MWD) are natural practical applications of the present invention. In a preferred 

embodiment of the invention, two accelerometers 3 and 4, are placed opposite each other, 
across the tool axis of rotation 5, as shown in Fig. 2 A. Accelerometer 3 provides 
orthogonal radial and tangential acceleration readings a^^, a^j and, optionally, the axial 
acceleration tool component a^i (not shown), Accelerometer 4 provides the radial and 

10 tangential acceleration readings a^, a^j, respectively, and optionally acceleration component 
a^j. Figure 2 shows an alternative preferred embodiment with foxir separate accelerometers, 
5, 6, 7, and 8, each providing a single acceleration reading, which together operate in the 
same manner as the aforementioned two accelerometers. The two accelerometers, 3 and 4, 
provide a^^, a^j, a^, and a^j, all located in the same plane. This arrangement provides a 

15 means of measuring lateral tool acceleration in the tool's rotating frame free of centriftigal 
and angular acceleration effects. In accordance with the present invention, a first and 
second magnetometer (exact position not shown) provide magnetic readings B^^ and By. The 
magnetometers are moxmted on or in the drilling tool 2 to provide the tool's orthogonal 
magnetic readings B^ and By relative to the earth's magnetic field vector. Preferably, the 

20 magnetometers should be placed on or in the tool linearly aligned with the position of the 
accelerometers and the tool axis of rotation 16. 

Additional parts of the tool of the drilling system are depicted in Figure 1 7, 
including the drill collar 17 and stabilizer 18. 

25 

B, The Method of the Invention 

In accordance with the preferred embodiment, the method of the invention 
pennits two independent corrections. The first correctioii removes the gravitational 
component from the acceleration readings that results when the tool is tilted away from the 
30 vertical direction. The second correction provides the lateral velocity of the drilling tool 
relative to a borehole reference frame. 

In a preferred embodiment, the method which corrects both inaccuracies 
includes the foUowing'steps (discussed again in greater detail later): 

35 
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(a) reading the instantaneous tool lateral acceleration components, a^,, a^z, a^,, and a^j, 
employing an accelerometer detection system and reading the tool's instantaneous magnetic 
values and By; 

(b) compensating for the centrifugal and radial acceleration components by 
5 calculating lateral tool accelerations in the rotating frame of reference: 

= (a^^ - ) / ^ ^^^^ acceleration in the direction of a,,; 

= - a^^)/ 2 the tool acceleration in the direction of a^^; 

(c) calculating the tooFs instantaneous magnetic phase relative to the earth's 
10 magnetic field (tool magnetic phase) (p^; 

(d) using the accelerometer and magnetometer instantaneous data readings, 
determining the phase shift cp^, which is the difference between the magnetic phase cp^ and 
gravitational tool phase (p, determining the borehole inclination gravitational component 
GsinCttj) relative to vertical; and 

15 (e) calculating the tool lateral acceleration components and optionally, correcting for 

the gravitational component, or converting the measurements to the borehole frame of 

reference, or doing both; and 

(f) calculating the initial velocity and the instantaneous velocity by integrating the 

acceleration components calculated in step (e). 
20 accordance with the present invention, the signals recorded by the 

accelerometers are related to other system variables by the following expressions: 

^ri = cos(fl? ) + siii(p ) + rfi) ^ + G sin(a .) cos(^ ) 
a ^2 = "^jc cos(^)-a^ sin(^ ) + r^y ^ - Gsin(a.)cos(^) 

dco 

25 = cos(^)+ a sin(§?)+ r— -+ Gsin(a Jsin(^) (1) 

at 



~dt 



a^2 = " cos(9? )- sm{(p ) + r —j^ + G sin(a ) sm(^ ) 



30 



35 



where: 

a^, ay, and a^ are the acceleration components of the tool's center of gravity relative 
to the Ugrehole XYZ reference fi:ame; 

(p is the instantaneous phase of the rotating tool ((p=0 when a^j is aligned with the X 

axis); 

ttj is the tool inclination angle in relation to the earth's gravity vector (vertical); 
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r is the rotational radius of the accelerometer; and 

G is the acceleration constant of the earth's gravitational field (s9.81 m/s^). 
Gsin(ai)sin((p) and GsinCOj )cos((p ) are the gravitational components arising from 
tool tilt away from vertical. 
5 The tool phase <p is: 

t 

9>(0 = 9>,o+ jo)it)dt (2) 

0 



where co is the instantaneous angular speed of the tool. From previous equations for a^i and 
aj2 the module of co is calculated as 



kl=J^ (3) 
J 2 and the angular acceleration is 

By tracking both \o) \ and dco/dt, co can be determined. 

20 Employing the above relationships, aiid the described hardware system and 

programming software, the method for obtaining lateral tool velocity with correction for the 
gravitational component and conversion of the velocity relative to a borehole reference 
frame is now discussed in detail. The method includes the following particular steps. 



25 Step a. Obtain instantaneous tool lateral acceleration components a^j, a^j, a^j, 

and a^, and the tool's instantaneous magnetic values and By. 

In a preferred embodiment, this step requires reading real-time data signals 
from the two (or more) accelerometers and two magnetometers to obtain the parameters a^j, 
30 3^2* at2, Bjj and By Parameters B^ and By are orthogonal magnetic phase readings relative 
to the earth's magnetic field. 



Step b. Compensate for the centrifugal and radial acceleration components by 

calculating lateral tool accelerations in the rotating frame of reference 
35 employing the following formulas: 
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= (a^i - ) / 2 the tool acceleration in the direction of a^,; 
= {a^^ - a^2)/2 the tool acceleration in the direction of a^,; 

5 Step c. Determine (p^,, the tool's instantaneous magnetic phase relative to the 

earth' s magnetic field (tool magnetic phase). 

The magnetic phase readings are used to determine the tool's magnetic phase 
with respect to the earth's gravitational pull. The direction of the magnetic field in space 

10 however, does not directly coincide with the gravitational pull; there is a phase difference 
(phase shift) of (p^. In most conditions, where the magnetic field disturbance is not strong 
and the borehole has a relatively constant direction, the phase shift (p^ will be a constant 
within the time firame of the few seconds necessary to determine the tool velocity. 
Therefore, the relationship (p = (pn, + (po where (p,, = a constant, can be reasonably assumed. 

15 Knowing B^^ and By, the tool's magnetic rotation phase (p^, may be obtained according to the 
present invention using the expressions: 

= 5sin(a„)cos(^„) 

5^ = -Bsin(a^)sin(^J 

20 

where B is the amplitude of the magnetic induction signal, and is the angle between the 
tool's axis and the earth's magnetic field vector. 

The tool magnetic phase <p„ is determined directly firom (5) provided that the 
borehole direction does not coincide with the direction of the B vector such that the noise 
25 level of the magnetic measurements is comparable to the signals B^ and By. Knowing B^ 
and By, the tool's magnetic rotation phase cp^j may be obtained by using the function 
<p = atan2(By3x) common to most mathematical function Ubraries. The function atan2 
resolves all four quadrants of the fiill angle (360 degrees). 

30 Step d. Determine the phase shift (p^ and the borehole inclination 

gravitational component Gsin(ai) relative to vertical, using the 
accelerometer and magnetometer instantaneous data readings, and 
calculate <p. 

35 
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If a correction for tool tilt is not desired, then it is unnecessary to determine 
Gsin(ai) in this step. However, it is the usual case to correct for the effect of tool tilt. The 
following procedure is used in a preferred embodiment to determine GsinCOj) and (p^, where 
G is the acceleration constant of earth's gravitational field (a 9.81 m/s^). The tool magnetic 
5 phase (p„ is known from the previous step. Gsin(ai) can be calculated uiider the assumption 
that the gravitational component does not contribute to the lateral acceleration of the tool. 

As shown in Figure 4, the acceleration signals, with centrifugal component 
removed, along with the (p„ phase information, are filtered through identical low-pass filters 
^ havi^g in a preferred embodiment a 20 ELz cutoff frequency. The 20 Hz cutoff is believed 
10 adequate to pass all gravity-related components, although it will be appreciated that different 
frequencies may be used in alternative embodiments. 

The signals are then decimated and fed into a quadrature detector known to 
those skilled in the art. In the quadrature detector both acceleration signals a, and a^ are 
multipUed by the sin((p^ and cos(<p^. The outputs averaged over time (few seconds in a 
15 preferred embodiment) yield two complex numbers c and d, where: 

/=1 -^^ 1=1. 



'imag 



25 



Where: N is the nimiber of signal samples processed during the averaging; 

a^ and a^j are consecutive samples of a^ and a^, respectively; and 

(pmi are consecutive samples of (p^j. 

Both complex numbers are 90 degrees out of phase since the gravitational 
component is 90 degrees out of phase in a^ and a^ respectively. The magnitude of these 
complex mmabers equals to 0.5 Gsin(ai) and the phase of c equals to (p,,, therefore: 



30 G sin(a,. ) = 2^c^^^ + 

(p^ = atan2(c,„,^,c,,,/) 

Once the phase shift (p^ is foimd from step (d), combined with the parameter 
(Pm known from the previous step, (p may be calculated according to the relationship: 
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The same information can be obtained from the complex number d, 
remembering that there is a 90 degree phase shift between c and d. If the magnitude and 
phase is obtained from both complex outputs, in a preferred embodiment it can be averaged 
to decrease uncertainty. 

5 This process yields both the phase shift (po and magnitude of the gravitational 

component Gsin(ai). The time constants of the averaging process can be as long as 30 
seconds or more, if the phase information from magnetic sensors is used, since there is no 
systematic drift between the (p^^ and (p^ other than changes of the borehole direction or of the 
magnetic field, which typically are very slow. 

10 To assess the quality of the real-time data, the standard deviation of each 

measiired/calculated quantity may be determined, if possible. If the same information is 
available from several sources, preferably the one with the lowest standard deviation is 
chosen. Based on individual uncertainty estimates, the uncertainty of velocity determination 
can be calculated and made available to the computer system for storage, 

1 5 While phase detection is desirably obtained by using magnetometers, this 

method is not available when the tool axis coincides with the magnetic vector. An 
alternative, although less accurate method of phase determination using the accelerometer 
signals, is available in accordance with a specific embodiment of the present invention. 
According to Eq. (2), the gravitational tool phase (p can be calculated as an integral of the 

20 instantaneous angular velocity co, which can be determined from Eq. (3) and Eq. (4). It will 
be appreciated that this approach is sensitive to accelerometer scale error and may suffer 
from poor resolution of© at low speeds. Nonetheless, in accordance with the invention, the 
approach can serve as a backup algorithm in situations where magnetic information is not 
available. 

25 

Step e. Calculate the lateral tool acceleration components in the borehole reference 
frame and, optionally, correct for the gravitational component, or convert the 
measurements to the borehole reference frame, or do both. 

30 To obtain lateral accelerations a^ and ay^ the raw acceleration signals are 

subtracted so that centrifiigal and angular acceleration components cancel out: 

^''^ J^''^ = cos(^) - ay sm(q>) + Gsm(a,.)cos(^) 
a ^ 

35 = sm((p)+ cos(^) + Gsm(a.) sm((p) 
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The signals above also contain the modulated gravitational component 
Gsin(ai)cos((p), Since Gsin(ai) and <p have been detemiined in the previous step, in 
accordance with the present invention the gravitational component can be subtracted from 
both signals yielding accelerations corrected for gravitational components a^ and a^g: 

5 cos(9?) - a sin(^p) = a 

\ (8) 
sin(^) + a cos(^) = 



10 



XYZ. 



Step f : Transforming the accelerations into the formation reference frame 

Solving Eq. (7) , yields: 

= cos(^) + a^^ sin(^>) 
a^^-a^ sin(??) + a,^ cos(9?) 

15 Equation (9) is used to convert the tool acceleration from the (r-t-a) reference 

frame to the XYZ borehole reference frame. All variables have been previously determined 
in order to calculate a^ and a^ Note also that Eq. (9) may be used when no correction is 
desired for the gravity effect of tool tilt on the accelerometers, and only a conversion to the 
borehole frame of reference is desired. 

20 

Step g. Calculate the lateral velocity components by calculating initial 

velocity and integrating the instantaneoxis acceleration found in step e. 

Knowing a^ and ay from the previous step, the lateral velocity components v,^ 
25 and Vy may be calculated. The lateral velocity calculation is provided in a preferred 
embodiment as follows: 



30 



35 



-13- 



NY2 - 1207701.1 



where v^^ and v^^y are uxiknown initial velocities at arbitrarily chosen time Tp. Since the 
borehole restrains the motion of the tool during any period, the lateral displacement is less 
than or equal to the slack As between the drill collar and the borehole wall. 

Since values of a^ and ay are known at any point in time, the initial velocities 
Vox 3nd Vpy can be calculated from: 



To 

V, 



'Ox 



10 



(11) 



15 



20 



with the uncertainty of the measurement method less than As/(d-To). For example, to 
achieve an uncertainty of 0.02 m/s in a borehole having a slack of 5 cm, the mioimum 
integrating time should be 2.5 seconds. 

After the individual lateral velocity components are extracted, the modulus of 
the lateral velocity is calculated in accordance with the invention as: 



25 



30 



35 



+ v'^ (12) 



Once the instantaneous velocity is calculated, the decision can be made whether to initiate 
an NMR measurement. 

In order to use the velocity calculation as described by equations (10-12) with 
computer processing, it is desirable to simplify the data processing to minimize the 
calculations. Thus, assxmiing a minimum of 2.5 seconds and a sampling frequency of 8 
kHz, the number of samples integrated would exceed 20,000. The memory requirement for 
direct implementation would be substantial. Therefore, in a preferred embodiment, a 
multiple-window approach is performed, wherein the integrals are calculated over K 
partially overlapping time windows as shown in Figure 5. The individual samples do not 
have to be stored, only the ^ptegrals and number of samples integrated. When an integrator 
reaches the preset nimiber of samples, i.e., 2.5 seconds worth of data in a specific 
embodiment, it becomes the source of velocity information for the system, until the next-in- 
line integrator reaches the minimvun number of samples. Then the first integrator is reset 
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and begins another new integration, while the second integrator provides velocity 
information. This processing approach tolerates some discontinuity in the velocity signal 
that is introduced when switching integrators in the Kth increase during processing. 
However, as simplified using the above approach the calculations are manageable and 
5 provide reasonably accurate results. The performance of recursive filters dxiring velocity 
. retrieval may also be tested in a specific embodiment. 

C. Tvpical Acceleration Signals Provided in the Tool Reference Frame and Svstem 
Bandwith Concems 

10 A data set obtained firom J. Dudley*s, *^MRIL MWD Sensor Physics" 

DNMWD016 Rev.la, Security DBS, was analyzed to assess the magnitude and spectral 
composition of typical lateral tool accelerations. This reference provides graphs of a 
complete data set in the reference firame of the rotating tool (drill bit), and is hereby 
incorporated by reference. The following example data analysis focuses only on the radial 

15 component of the tool acceleration and only for selected firagments. No caUbration 
information is provided. 

The centrifiigal acceleration depends on the radius and angular velocity of the 
tool. The accelerometer's operating radius is limited by the tool's diameter. The average 
angular velocity of the tool equals the driving velocity. However, under stick-slip 

20 situations, the instantaneous velocity may change quite rapidly. 

The natural frequency of the tool string and forces on the bit modulate the 
instantaneous angular velocity of the tool. Even during relatively vibration-firee periods, o 
tends to oscillate, as illustrated in Figure 6. The angular velocity modulation is much more 
dramatic in stick-slip conditions, where the drill bit actually stops rotating for firactions of a 

25 second. From a momentary standstill, the bit accelerates when the moment exerted by the 
winding drill pipe exceeds the static fiiction. The peak angular velocity under stick-slip 
conditions is more than twice the average velocity. Figure 7 illustrates the condition. 

The acceleration data provided by Sperry-Sim in a report by J.D. Jansen, 
"Whirl and Chaotic Motion of StabiUzed Drill Collars" SPE 20930, pp. 435-448, are in the 

30 rotating tool reference frame. The reference is hereby incorporated by reference. Figure 8 
illustrates the Fast Fourier Transform (FFT) of two 16 second fragments of acceleration 
data. The spectrum contains strong peaks of frequency modulated base rotational 
component caused by the effect of gravity on the tilted tool. The peaks are spaced about 0.7 
Hz apart, which is the modulating frequency likely closest to the natural frequency of the 

35 tested toolstring. For comparison, a signal was modeled with 0.92 Hz base frequency 
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modulated with 0.67 Hz non-harmonic signal resembling the measured signal, as illustrated 
in Figure 9. Its FFT is shown in Figure 10. The modeled signal has a spectrum resembling 
that of the measured signal. However, the ampUtudes of the individual peaks, which are 
govemed by Bessel functions, depend on a deviation coefficient. Both FFT's in this 

5 example were calculated from a set of 4096 signal data points sampled at a rate of 250 Hz. 

The frequency spectrum in Figure 7 displays both the gravitational 
component dominating at lower frequencies, and the lateral motion component. An ideal 
signal would display only the gravitational component. The gravitational component is 
bandwidth limited to approximately 2(Af+24,), where Af is the maximum frequency 

10 deviation and 4i is the maximum modulating frequency. Assuming a S Hz maximum 
frequency deviation during drilling and a maximum modulating component of 2.5 Hz, the 
bandwidth of the gravitational signal would stretch up to 20 Hz. 

Since the lateral motion of the tool can be described in the frequency domain 
as a sum of harmonic motions of various frequencies, it is possible to determine the . low 

15 frequency cutoff of the measurement system based on a desired velocity accuracy, and a 
maximimi allowable amplitude of the tool's movement which is limited by the borehole 
walls. See Table 1 . For example, for a typical 2 cm or 4 cm peak-peak amplitude 
limitation, the system must pass accelerations down to 0.2 Hz in order to achieve a precision 
of 0.025 m/s. The peak acceleration in that case would be approximately 3 mG, this value 

20 giving an indication of the required maximvim resolution of the accelerometer at the lowest 
frequencies. 



Table 1. Max velocity vs. frequency and amplitude of harmonic lateral motion 



25 



30 



'4 

35 



AinpBtiide [mp» 


0.0001 


0.0002 


0.0005 


0.001 


0.002 


0.005 


0.01 . 




0.05 


0.1 


f 


m 
























[rad/,1 










Nfax velocity in m/s 










0.1 


0.63 


O.OOOl 


0.0001 


0.0003 


0.0006 


0.0013 


0.003 


0.006 




0.031 


0.063 


0.2 


1.26 


0.0001 


0.0003 


0.0006 


0.0013 


0.0025 


0.006 


0.013 




0.063 


0.126 


0.5 


3.14 


0.0003 


0.0006 


0.0016 


0.0031 


0.0063 


0.016 


0.031 




0.157 


0.314 


1 


6.28 


0.0006 


0.0013 


0.0031 


0.0063 


0.0126 


0.031 


0.063 




0.314 


0.628 


2 


12.57 


0.0013 


0.0025 


0.0063 


0.0126 


0.0251 


0.063 


0.126 




0.628 


1.257 


5 


31.42 


0.003 


0.006 


0.016 


0.031 


0.063 


0.157 


0.314 




1.571 


3.142 


10 


62.83 


0.006 


0.013 


0.031 


0.063 


0.126 


0.314 


0.628 




. 3.142 


6.283 


20 


125.66 


0.013 


0.025 


0.063 


0.126 


0.251 


0.628 


1.257 




6.283 


12.566 


SO 


314.16 


0.031 


0.063 


0.157 


0.314 


0.628 


1.571 


3.142 




15.708 


31.416 


too 


628.32 


0.63 


0.126 


0.314 


0.628 


1.257 


3.142 


6.283 




31.416 


62.832 


200 


1256.64 


0.126 


0.251 


0.628 


1.257 


2.513 


6.283 


12.566 




62.832 


125.664 




314159 


0.314 


- 


1.571 




<5 283 


15.708 







157.080 
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1000 



I 628319 I 0.628 \ 1.257 | 3.142 | 



6.283 













12.566 


31.416 


62.832 




314.159 





0.0001 1 


0.0002 1 


0.0005 1 


0.001 1 


0.002 1 


0.005 1 


0.01 1 


0.02 


0.05 1 O.l 


f 
























rHzi 


trad/*} 










Max Acceleration in G 










0.1 


0.63 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.001 


0.002 


0.004 


0.2 


1.26 


0.000 


0.000 


0.000 


0.000 


0.000 


0.001 


0.002 


0.003 


0.008 


0.016 


0.5 


3.14 


0.000 


0.000 


0.001 


0.001 


0.002 


0.005 


0.010 


0.020 


0.050 


0.101 


1 


6.28 


0.000 


0.001 


0.002 


0.004 


0.008 


0.020 


0.040 


0.080 


0.201 


0.402 


2 


12,57 


0.002 


0.003 


0.008 


0.016 


0.032 


0.080 


0.161 


0.322 


0.805 


1.610 


5 


31.42 


0.010 


0.020 


0.050 


O.lOl 


0.201 


0.503 


1.006 


2.012 


5-0 


10.1 


10 


62.83 


0.040 


0.080 


0.201 


0.402 


0.805 


2.012 


4.024 


8.049 


20.1 


40.2 


20 


125.66 


0.161 


0.322 


0.805 


1.610 


3.219 


8.049 


16.1 


32.2 


80.5 


161.0 


50 


314.16 


1.01 


2.01 


5.03 


10.06 


20.12 


50.30 


100.6 


201.2 


503.0 


1006.1 


100 


628.32 


4.02 


8.05 


20.12 


40.24 


80.49 


201.2 


402.4 


804.9 


2012.2 




200 


1256.64 


16.1 


32.2 


80.5 


161.0 


321.9 


804.9 


1609.7 


3219.4 






500 


3141.59 


100.6 


201.2 


503.0 


1006.1 


201^2 


5030.4 














402.4 


804.9 


2012.2 


4024.3 


8048.6 













The above estimate is only valid for a system that does not rotate. If the tool 
rotates with constant cd^, then the acceleration induced by the slow lateral motion co would 
be modulated as the rotational signal shifts to (o^-co and co^,+co. However, there are 
conditions under which the tool stops rotating while performing measurements. 

Table 2 above shows the peak accelerations observed during harmonic lateral 
motion. For example, to observe 50 G acceleration, the tool must vibrate at a frequency of 
50 Hz with an amplitude of 5 mm or 10 mm peak-peak. 

The strongest accelerations were observed during bit whirl. Peak values of 
about 13-14 G were detected in Dudley's sample accelerations. As shown in Figure 11, the 
FFT of the accelerometer signal during bit whirl indicates high frequency components 
peaking around 15 Hz with substantial content extending past 100 Hz. The frequency 
beyond 100 Hz in this sample is xmreliable owing to^ a 250 Hz sampling rate and additional 
use of a 125 Hz low pass filter for internal filtering. If resolution of higher frequencies in 
excess of 100 Hz is desired, data should be collected at faster sampling rates without digital 
filtering. Additionally, it is advisable to shape the frequency response of the anti-aliasing 
filter to prevent fiirther distortion of the spectrum. It should be noted that the observed peak 
acceleration will be affected by the spectral response of the entire measuring system. 
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Close et al., lADC/SPE 17273, pp. 559-568, conducted measxirements of 
Borehole Assembly OBHA) vibration while drilling using a sampling rate of 2 kHz. The 
reference is incorporated herein by reference. The analysis revealed peak lateral 
accelerations of 13 G during reaming, 2G behind the positive displacement downhole 

5 motor, and 25 G while drilling through casing shoe. Accelerations in excess of 100 G were 
also reported when the BHA rotational speed matched the BHA's lateral resonant frequency. 
Accelerations of that magnitude are frequently destructive. 

From published data it is apparent that shocks in excess of 50 G are rare and 
are not considered within normal operating conditions. Therefore, the accuracy of a motion 

10 detection system under these infrequent conditions is not critical. However, the system 
should be able to recover from such severe shocks without lasting adverse effects. 

To estimate the acceleration magnitude and shape during collision of the tool 
with a borehole, a simple model was used, simulating a section of the drill collar colhding 
with a perfectly rigid barrier. The elasticity K was approximated based on the drill collar 

15 size and the material property. The simulated acceleration at the pipe's center is shown in 
Figure 12. The pulse duration is independent of the initial velocity and is affected only by 
the pipe's properties. As shown in Figure 13, the peak acceleration ampUtude increases 
Unearly with the initial velocity. This simulation of the collision can be considered a worst 
case since, in reality, the borehole fluid and cmshing of the rock cushions the impact. 

20 Nonetheless, the simulated signal provides a good base to calculate the effect of bandwidth 
limitations on the accuracy of velocity. 



Table 3. Velocity and 0.5 ms displaciement errors vs. bandwidth (50 G maximum 
acceleration, simulated elastic collision) 



25 



bandwidth (-3Db) 


velocity error 


0.5 ms displacement error 


FHzl 


[m/sl 


[mm] 


200 


0.604 


0.297 


500 


0.369 


0.18 


1000 


0.188 


0.09 


2000 


0.093 


0.045 


5000 


0.036 


0.018 



Figure 14 shows the acceleration signal after passing through a 500 Hz low- 
35 pass filter. The group delay of the filter causes a time shift between the input and output 
signals. Consequently, the calculated velocity will be loaded with an error as shown in 
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Figure 15. No post-processing can compensate for that group delay error and an accurate 
tool velocity in real-time cannot be obtained. Therefore, an adequate bandwidth for the 
accelerometers, as well as for the input anti-aUasing filters are chosen in advance in 
accordance with a preferred embodiment. If, on the other hand, a correction of the NMR 

5 signal based on actual velocity is considered, then the group delay can be digitally 
compensated employing phase correctors. 

Estimates provided in Table 3 and Figure 16 indicate that a bandwidth of 4 
kHz is required to assure a velocity measurement precision of 0.05 m/s during impacts of 50 
G or, equivalently, as found in simulation, a 0.35 m/s impact velocity. The velocity error is 

10 proportional to peak acceleration during impact, 

D. Uncertaintv and Error Analvsis 

The accuracy of velocity calculations based on acceleration measurements 
are affected by both instrumental factors and motion characteristics. Some factors produce 

15 purely systematic errors, such as gain error, while others, despite their systematic character, 
are exhibited as semi-random errors owing to a randomizing effect of the conditions, for 
example, with bandwidth effect and quantization noise. Some factors, such as transducer 
noise are purely random. Owing to their dual nature, some uncertainty components may be 
described both statistically and as a worst case. 

20 The effect of quantization noise can be approximated by a random noise with 

a variance of l/12*q^, where q is the least square best fit value with a conversion resolution 
in m/s^. Since v^^ and Vy are integrated over the period T containing T/dT samples, with dT 
the sampling interval, the variance of the velocity is: 

< = < = {[h"^] If = [-k"') ^""^ 

The factor of V2 is a result of averaging the outputs of two accelerometers to obtain 
compensated acceleration signal (1 1). For a 14 bit+sign A/D converter and 50 G range, q is 
0.03 m/s^. Assuming a T of 2.5 seconds and a dT of 125|is, the variance of the individual 
30 velocity components is 1.17*10 ^ resulting in a standard deviation of 0.0001 1 m/s. To 
simplify the calculations, the variance of total velocity v, as described by Eq. (14), can be 
estimated to be less than 2a^V^. The variance is proportional to velocity integration time 
and sampling interval and to the square of the A/D converter's resolution. Each additional 
bit reduces the standard deviation by a factor of 2. 

35 
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The transducer noise Oj can be transformed the same way as the quantization 
noise according to the expression: 




5 

except the input level is bandwidth dependent- Assuming a noise level of 1 mG/sqrt(Hz) 
using a Silicon Designs accelerometer having a bandwidth of 4 kHz, the input noise level 
can be estimated at 0.62 m/s^ rms. Applying Eq. (15) and multiplying the result by 2 yields 
a velocity variance of 0.00012 mVs^ and a standard deviation of 0.011 m/s. 

10 Any calibration error for the accelerometers will have a proportional effect 

on the calculated velocity. Since, over a long period of time, the average acceleration is 
zero, the effect of gain error is not cumulative. For any time period this error is proportional 
to the average acceleration for that time and the time duration. As a first approximation, the 
direct coUision with the borehole wall may be analyzed. For the duration of the collision, 

15 the average acceleration is approximately 30 G and the pxilse duration is 2.2 ms. The 
velocity error accumulated during the collision resulting from a 1% calibration error is 
0.0065 m/s. The gain of the accelerometers can be initially calibrated to a measurement 
precision better than 0.5% using the earth's gravitational field as a reference. 



20 temperature coefficient of gain in the order of 2% over a 125 °C temperature span. This 
temperature behavior is repeatable so that compensation is possible. It is realistic to assume 
a net temperature coefficient after compensation to be less ihzn 0.5%. Similarly, the 
resulting velocity error would be 0.0032 m/s. 



25 freedom the tool has in the borehole. The larger the borehole compared to the tool's 
diameter, the larger the error according to Eq. (1 1). 



Commercially iavailable micro-machined siUcone accelerometers have 



The limited integration time results in an error dependent upon the degree of 



Any fluctuation of phase information used to compensate the gravitational 
field can cause a disturbance in the accelerometer information since the gravitational 
component is removed by subtracting the following term froin equation (7): 



30 




The sensitivity of this correction factor to phase noise is: 




corr 



dq> 



= Gsin(a,.)sin(^) 



(17) 



35 
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which peaks at Gsin(ai). Consequently, the variance of the acceleration signal due to the 
phase noise is: 

gI « G^sm\a,)Gl (18) 



10 



15 



Similarly as in IJq. (14) the resulting velocity variance due to phase noise is: 



(19) 



The noise approaches zero when approaches zero, representing a perfectly vertical 
borehole. 

Fluctuations of the phase signal come from magnetometer noise, extemal 
magnetic field disturbances and fluctuations of the tool's rotational axis in relation to the 
borehole axis. System problems occur from disturbances and fluctuations of the tool's axis 
in relation to the borehole axis. The use of magnetometers to detect phase is problematic 
when the direction of the borehole approaches the direction of the earth's magnetic field 
vector. 

Since the phase is calculated as: 



20 



= -arctan^^j. 



1 



1 



- — < — ^ < 
2-5,-2 



(20) 



25 



the sensitivity of O^, to a B disturbance is: 



SI:- 



1 



1 




(21) 



30 



which assumes its lowest value of when By =0. At this point, the value of B^ is sin(0, 
where a„ is the angle between the borehole and the earth's magnetic field vector. A 
symmetrical relationship exists for B, noise, but a reasonable estimate can be based on just 
one magnetic field component since the influence of one component peaks when the 
influence of the other is the lQ;^est. The phase noise caused by the magnetic signal noise 
can be expressed as: 



35 
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^2 _ ^2 
^9 = ^9. 



(22) 



According to Eq. (22), the phase noise increases as goes to 0, representing drilling along 
the earth's magnetic field vector. Combining Eq.(19) and Eq.(22) a formula linking 

» - 

magnetic field noise with velocity is derived: 



10 



_ ^2 _ J: 



T G^sm^ia,)TdT 



(23) 

The noise of the magnetic induction signal has three uncorrelated components: (1) 
magnetometer noise, which is a function of the device used to measure the magnetic field; 
(2) fluctuations of the tool's axis in relation to the earth's magnetic field; and (3) influence 
of the antenna-magnet. The noise of the magnetoresistive sensor, in this case a Honeywell 
HMC1002, for a bandwidth of 0.01-20 Hz, is under 6 uGauss rms with a variation of 
36*10^2 Gauss^ 

The tool's axis is constrained by the borehole and, assuming a 30 fl: distance 
between stabilizers and maximum 2 inch slack between the borehole wall and the stabilizer, 
as shown in Figure 17, the maximum ampUtude of p is 2*arctan(2"/30')=0.64°. Assuming a 
uniform distribution of P within ±p/2 range, the variance of this distribution is pVl2. For 
small P the variance of By due to the tool's axis fluctuations is: 



20 



25 



30 



35 



(24) 



where the variance of P is in radians^. 

The magnetic field generated by the NMR magnet may produce additional 
noise if the position of the magnet changes in relation to the magnetometer owing to the 
deformation of the tool's structure or the magnetometer's mounting. Only components 
within the 20 Hz bandwidth are significant since higher firequencies will be filtered out. 

Figure 18 illustrates the effect of the antenna magnet on the magnetic sensor 
measurements. When the magnetometer is on the magnet's axis it registers some By field 
induced by the magnet. Magnet-induced components in other directions are zero in that 
configuration. Figure 19 shows the relationship between the magnet-magnetometer distance 
and the By component. At a distance of 70 inches, the magnet's field equzds the earth's 
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magnetic field. The distance L should be chosen such that the magnet's field, together with 
the earth's field, do not saturate the magnetometer at any moment. A distance of 100 inches 
appears adequate in the specific embodiment. 

When the tool is bent by an angle y, the By component will change slightly. 
5 Figure 20 shows the sensitivity of the By to y for up to a 1 deformation. Assuming a 
standard deviation of y not exceeding 1 the resulting By variance is 1.9*10"^ Gauss^. For 
the drill collar, the maximum specified bending radius is 18 ''/1 00 ft. At a distance of 100 
inches the sensitivity is about 0.0025 Gauss/radian. 

10 

Table 4. Random error budget 



Source 


Description 


Uncertainty contribution 
[la/s la] 


Quantization noise 


14 bit+sign A/D converter, ± 50G range, 8kHz 
sampling rate, 2.5s averaging time 


0.0002: 


Transducer noise 


4kHz bandwidth and 6mG/sqrt(Hz) noise level 


0.011 


Phase detection error 
(antenna magnet effect) 


Magnetometer 100" from magnet's center, 
a„>1.5- 


0.03 




Total (quadrature sum) 


0.032 



20 

Table 5. Random error budget 



25 



30 



Source 


Description 


Uncertainty 
contribution [m/s 
lol 


Gain calibration error 


A 50G peak acceleration during a wall collision 
2.2ms long; 0,5% calibration error 


0.0032 


Gaia temperature effect 


As above, assuming 0.5% net gain change after 
condensation over the working tenperature change 


0.0032 


High frequency cutoff 


4000Hz bandwidth, tnavi'mnm 50G acceleration (see 
Table 3) 


0.05 


Limited integration time 


Assuming 5cm slack in the borehole and 2.5s 
integration time (see formula (14)) 


0.02 




Total (worst case) 


0.0764 



A combined effect of all 3 magnetometer-related effects calculated according 
to (23) is shown in Figure 21. The velocity uncertainty (la level) is below 0.015 m/s if the 
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direction of the borehole is at least 1.5° apart from the direction of the earth's magnetic field 
vector. Otherwise, if the borehole direction is essentially parallel to the earth's field vector, 
the phase information has to be derived from the acceleration signals themselves. 

Tables 4 and 5 simmiarize the uncertainty budget in accordance with a 
5 specific embodiment. It will be appreciated that the major contributor of random noise is 
the phase noise. The systematic error is dominated by the bandwidth limit. In accordance 
with the present invention, this error can be reduced in post processing by introducing phase 
correctors to reduce the group delay of the filters. 

10 It will be readily apparent to those skilled in the art that various 

modifications may be made without departing from the spirit and scope of the invention. 
The scope of the invention shall not be limited to the embodiments described in the 
specification but shall be defined by the claims as appended. 
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